Noise Reduction and Universality in Limited Mobility Models of Nonequilibrium 

Growth 
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We show that a multiple hit noise reduction technique involving the acceptance of only a fraction 
of the allowed atomistic deposition events could, by significantly suppressing the formation of high 
steps and deep grooves, greatly facilitate the identification of the universality class of limited mobil- 
ity discrete solid-on-solid conserved nonequilibrium models of epitaxial growth. In particular, the 
critical growth exponents of the discrete one dimensional molecular beam epitaxy growth model are 
definitively determined using the noise reduction technique, and the universality class is established 
to be that of the nonlinear continuum fourth order conserved epitaxial growth equation. 



PACS: 05.40.-l-j, 81.10.Aj, 81.15.Hi, 68.55.-a 

In 1991 Das Sarma and Tamborenea introduced j|] 
an extremely simple one dimensional (d= 1-1-1) instan- 
taneous relaxation limited mobility conserved discrete 
solid-on-solid model of ideal molecular beam epitaxial 
(MBE) growth under random vapor deposition nonequi- 
librium growth conditions. In spite of the decep- 
tively simple deposition and relaxation/incorporation 
rules controlling its growth dynamics the universality 
class of this discrete growth model, particularly in d=l-|-l 
dimensions, has remained controver sial Ig and unre- 
solved in spite of a substantial body of work |^-|ll[l . Our 
lack of understanding of the universality class of this 
one dimensional growth model |jl| is particularly mys- 
terious for the following three reasons: (1) Recent large 
scale simulations of the correspending two dimen- 
sional (d=2-|-l) growth model seem to fairly unambigu- 
ously indicate the 2-1-1 - dimensional growth universal- 
ity class to be that of the fourth-order nonlinear con- 
served MBE growth equation ||l^; (2) a number of the- 
oretical approaches based on the kinetic master equa- 
tion technique p"^-p^ as well as symmetry arguments |^] 
lead to the conclusion that the one dimensional model 
should belong to the fourth-order nonlinear conserved 
MBE growth equation; (3) extensive large scale simu- 
lations (using upto lO^** deposited atoms in the largest 
simulations) in d= 1-1-1 produce |^-^ excellent scaling of 
the dynamically evolving surface roughness with the scal- 
ing exponents, however, being approximately consistent 
with the linear fourth-order conserved growth equation 
p3| rather than the expected nonlinear one (with the 
additional complication of there being substantial 
skewness in the growth morphology implying that a lin- 
ear description cannot apply). In addition, the model 
exhibits an intriguing anomalous multiscaling |^-|l^ be- 
havior in the height correlation functions, which tran- 
scends the standard self-affine dynamic scaling ansatz. 
In this paper we obtain the elusive "correct" asymp- 
totic universality class of this one dimensional |l) discrete 
growth model by introducing a multiple hit noise reduc- 
tion scheme which has earlier been successful |l^,|l^ in 
the identification of the growth universality classes for 
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FIG. 1. (a) Dynamical morphologies from the system of 
substrate size L = 1000 at 10'*, 10^ and 10*^ monolayers (ML) 
for the original model, ie. m = 1. Inset: Schematic config- 
uration defining growth rules in 1-1-1 dimensions, m hits are 
needed at a site for an actual deposition event (m = 1 in 
ref. 1). (b) Morphologies from the m = 10 model plotting 
on the exact same scale as in (a) showing the much smoother 
surfaces. The small insets show the morphologies at 10'* and 
10® ML in appropriated expanded scales so that the detailed 



rough morphology can be seen. 
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FIG. 2. The interface width W{t) as a function of depo- 
sition time. Solid lines indicate the power law fit with the 
growth exponent (3 for system L = 1000 with (a) m = 1 and 
(b) m = 10. Inset: Plots of logj^Q Wa vs logj^Q L where Ws is 
the saturation width. Slopes yield the roughness exponent a. 
Right inset in Fig. 2 (a): W{t) for WV model with m = 1, 
5, and 15 from top to bottom. In the original model (m — 1) 
/3 ~ 0.36 whereas in the noise reduced (m = 5, 15) results 
13 ^ 0.26. 



Eden and ballistic deposition models. Our noise reduc- 
tion results establish that the discrete limited mobility 
growth model introduced in ref.l does indeed belong to 
the fourth-order nonlinear conserved MBE growth equa- 
tion. 

Our growth model |Q] is shown in Fig. 1, where we also 
show the dynamical growth morphologies for the original 
model and the noise reduced model. The noise reduc- 
tion scheme rescales time by the noise reduction factor 
m (where m is the number of attemps required at a site 
for an actual deposition process to occur — m = 1 in 
the original model of ref.l), and all our times are given 
in terms of this rescaled time (which also defines the av- 
erage film thickness in our model). The most striking 
aspect of Fig. 1 is that noise reduction is seen to have 
a drastic effect on the growth morphology — it strongly 
reduces the high surface steps and deep grooves which 
were the hallmark of the original limited mobility 
growth model. This suppression of high steps (equiva- 



lently, the reduction of local slopes) in the growing sur- 
face is, in fact, the key to the success of the noise reduc- 
tion scheme in obtaining the asymptotically correct uni- 
versality class of the growth model. We emphasize, how- 
ever, that noise reduction, while drastically suppressing 
surface high steps/deep grooves, still produces a skewed 
growth morphology where the up-down symmetry of the 
starting flat substrate is spontaneously broken by the 
nonequilibrium growth process — in fact the measured 
skewness in the growth morphology is found to be in- 
dependent of the noise reduction scheme. The skewness 
(si —0.5) in the growth morphology (even in its saturated 
steady state) is a unique characteristic of the nonequilib- 
rium growth model of ref.l which is not shared by the 
other limited mobility epitaxial growth models ||l9| , po|] ex- 
isting in the literature. 

To proceed further we use the dynamical scaling 
ansatz, and discuss the evolving surface kinetic rough- 
ness in terms of two independent critical exponents /? 
and a. The interface width or the root mean square fluc- 
tuation in the dynamical surface height h{x,t), where h 
is the height of the growing surface at (reduced) time t 
at the substrate spatial point x, is defined as W{L,t) = 
{{h — {h})'^y^^ with the angular brackets representing 
average over the substrate of width L (in the one dimen- 
sional X direction) as well as a noise ensemble average. 
The dynamic scaling ansatz (which is obeyed extremely 
well by our results for all values of the noise reduction fac- 
tor m = 1 - 15 we investigated) asserts that W{L, t) ~ 
for ^ < L and W{L,t ^ oo) = Ws{L) ~ L" in the sat- 
urated steady state for ^ 3> i, where the lateral corre- 
lation length ^ = ^{t) ~ t^^^. The growth(roughness) 
exponent (s) /9(a) and the dynamical exponent z = a/P 
(which describes the approach to the steady state associ- 
ated with the growth of lateral correlations) completely 
define the universality class of the growth model provided 
the growth problem is self-affine. In Fig. 2 we show our 
calculated W{t) and Ws{L) dynamic scaling plots for var- 
ious values of the noise reduction factor. While all the 
results show excellent dynamic scaling, it is clear that 
the critical exponents for the noise reduced model are 
/3 « 0.33, a « 1.0, 2 « 3 whereas the original model 
(to — 1) gives /? « 0.37, a « 1.4, z w 3.9 in agreement 
with earlier findings p|,^[7|-pT|] . We believe that the crit- 
ical exponents f3 ~ 0.33, a « 1.0, z « 3 are the correct 
asymptotic exponents defining the universality class of 
the one dimensional growth model originally introduced 
in ref.l, and that the noise reduction scheme successfully 
suppresses the correction to scaling which dominates the 
TO = 1 version of the model for many decades in time. 

The coarse grained continuum equation which is be- 
lieved |]7|"|ll]| to describe the growth model of ref.l is the 
conserved nonlinear MBE growth equation [n3| : 
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where h = h{x,t) is now the height fluctuation, h — (h), 
around the average surface height, and 77(2;, t) is the 
(white) shot noise associated with the deposition beam 
fluctuations which produce the kinetic surface roughen- 
ing. The critical exponents for the corresponding fourth 
order linear equation |p],^, where A4 = = 0, are 
trivially known to he P — 0.375, a — 1.5, z = 4. It was 
already noted in ref.l that the simulated growth expo- 
nents of the discrete model seem to be following those of 
the linear version of Eq. (|l|) although no particular signif- 
icance was attached to this fact. The intriguing aspect 
of the critical behavior of the discrete model in d= 1-1-1 
dimensions has been its consistency with the linear ver- 
sion of Eq.(|l]) as far as the global exponents, /3 (~ 0.375), 
a (~ 1.5), and z (~ 4) go whereas at the same time 
the up-down symmetry of the growth problem, which is 
manifestly present in the linear equation because h —h 
leaves the equation invariant, is broken with the evolving 
growth morphology explicitly showing a finite skewness 
s = (/i3)(/i2)-3/2 ~ _o.5 (obviously s ee for the lin- 
ear equation). Thus the puzzle until our current work 
has been that the one dimensional discrete growth model 
should not belong to the fourth order conserved linear 
growth equation universality, except that it does for as 
long as (at least upto eight decades in time) one can dy- 
namically simulate the model [iHl^]. Our multiple hit 
noise reduction scheme resolves the mystery by obtain- 
ing the asymptotic exponents by successfully eliminating 
the severe correction to scaling problem which hinders 
the 171 = 1 version [Q of the model. This is similar to 
what was earlier found in the Eden model p^ , p^ . 

The critical exponents we obtain in the noise reduced 
model (/3 ~ 0.33, a ~ 1.0, z ~ 3) are consistent with the 
1-loop dynamical RG |]l3| treatment (/3 = 1/3, a = 1, 
z = 3) and direct numerical simulations 0,|l| of the 
fourth order nonlinear conserved MBE growth equation 
HI with \2n ^ 0, but 1/4, A4 7^ in Eq.(0). The origi- 
nal thinking jT^ that these one loop results may in fact 
be exact for the fourth order nonlinear equation has re- 
cently been questioned with a 2-loop dynamical RG 
treatment obtaining miniscule (less than 0.5%) numerical 
corrections. These 2-loop corrections are too small to be 
of any practical significance to our work (or to other sim- 
ulations and experiments). We therefore conclude that 
the critical exponents of the noise reduced version of the 
growth model introduced in ref.l are the same as those of 
Eq.(|l]) with = for n > 3, and therefore the model 
belongs to the universality class of the nonlinear fourth 
order conserved MBE growth equation. 

We believe that the success of noise reduction in pro- 
ducing the correct asymptotic growth universality arises 
from the suppression of the infinite series of nonlinear 
term (A2„ ^ 0) in Eq.(P which are associated with the 
high steps in the original (m = 1) discrete growth model 
[01 . The role of the infinite series of nonlinear terms (i.e. 
\2n ^ for n = 3,4, ...) in Eq.(|l|) is quite subtle @]ll) 
because the existence of such an infinite series of relevant 
terms is unusual in critical phenomena. High steps, how- 
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FIG. 3. (a) Anomalous multiscaling behavior of the height 
correlation functions Gq{r,t) = {\h(x -|- r,t) — h{x,t)\'')^^'' at 
fixed t = 10^ ML and the nearest-neighbor height difference 
correlation functions Gq{l,t) = {\h{x + l,t) - h(x,t)\'')^''^ 
(inset) with g = 1 — 4 from bottom to top. Substrate size 
L = 1000 and m = 1. (b) The same plots for system with 
m = 10. The noise reduced correlation functions show only 
very weak anomalous multiscaling behavior. 

ever, imply large values of the slope dh/dx, indicating 
the existence of the infinite series. Simple power counting 
shows all the terms in this infinite series to be marginally 
relevant operators in the one dimensional growth prob- 
lem, and since they are all allowed by the symmetry of 
the discrete model, they should all be present |0,|ll[ in 
the growth model of ref.l. Note that the corresponding 
linear problem has a (= 1.5) > 1, implying that the in- 
finite nonlinear series should be generated if allowed by 
symmetry. This infinite series has recently been shown 
pO| , pl| to be responsible for giving rise to the anomalous 
"intermittent" multiscaling behavior in the growth model 
of ref.l which has attracted considerable attention [p|-pT|. 
In Fig. 3 we show that the noise reduction technique es- 
sentially eliminates the anomalous multiscaling behavior 
in the height correlation functions of the growth model 
of ref.l. Thus the noise reduction technique resolves all 
three of the intriguing and puzzling features of the one di- 
mensional minimal MBE growth model of ref.l, namely, 
(1) it establishes beyond any reasonable doubt the correct 
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universality class of the model to be the 'expected' fourth 
order conserved nonlinear continuum growth equation; 
(2) it suppresses the "unphysical" feature of high steps 
and deep grooves in the growth morphology while main- 
taining a finite skewness in the morphology; (3) it sup- 
presses the anomalous multiscaling in the model. 

Before concluding we emphasize an important feature 
of our growth model which is characterized by what is 
absent from Eq.(|l|) — the well-known Lapla- 
cian i^2§5T term associated with the generic Edwards - 
Wilkinson (EW) universality class is strictly absent from 
the growth equation describing the discrete growth model 
of ref.l. The model of ref.l is, in fact, the only known 
limited mobility MBE growth model which does not be- 
long to the generic EW growth universality class. It is 
worth emphasizing this point because this has been a 
controversial and contentious [^,^,^ issue in the past. 

To reinforce this point we have also carried out noise 
reduction simulations of the closely related Wolf- Villain 
(WV) model which differs from the model of ref.l 
only in that all deposited atoms, independent of their ini- 
tial coordinations, are allowed to move to lateral nearest- 
neighbor sites in order to maximize their local coordina- 
tion number. Although it is well-accepted [pp^] that the 
WV model asymptotically belongs to the EW universal- 
ity class, this crossover has never been clearly observed 
in simulations because for all practical purposes the dy- 
namic scaling behavior of the WV model is similar 
to that of the model of ref.l upto the longest simula- 
tion times. Our d=l-|-l noise reduced WV model simu- 
lations, shown as an inset in Fig. 2 (a), clearly show that 
the asymptotic growth exponent (3 decreases from ~ 0.36 
(for m = 1) to ^ 0.26 ~ (3ew in d=l-|-l (for m = 15) un- 
der the noise reduction technique. Thus the WV model 
belongs to the EW universality and the model of ref.l 
belongs to the fourth order nonlinear MBE growth uni- 
versality. 

The non-existence of the EW term (which, if it existed, 
would have defined the universality class of the model be- 
cause all the fourth order terms in Eq.(|^) are irrelevant 
compared with the EW Laplacian term) in the growth 
model of ref.l is, in fact, an exact result due to a hid- 
den symmetry in the growth model which produces 
an exactly vanishing surface current [ p5[ in the model on 
a tilted substrate. The surface current on a tilted sub- 
strate is, in general, proportional [^,^ to the strength 
1^2 of the generic EW term in the growth equation, and 
therefore its vanishing implies 1^2 = 0. Our calculated in- 
clination dependent surface current on a tilted substrate 
in the growth model is essentially zero within error bars 
for all values of the noise reduction factor (m = 1 - 15). 
This is also a strong indication that noise reduction does 
not change the universality class of our growth model. 
This is similar to what was found earlier |l^,|l^ in the 
Eden model. It may be worthwhile to point out that in 
the Eden model also 0,11^ the success of noise reduc- 
tion in eliminating corrections to scaling arises from the 



suppression of high steps. 

We conclude by discussing why understanding the uni- 
versality class of the growth model introduced in ref. 1 is 
of considerable interest. An important theoretical rea- 
son is that this model is the only limited mobility MBE 
growth model which does not belong to the generic EW 
universality class, and therefore, as an exception, its 
proper theoretical understanding is of obvious interest. 
The fact that this growth model exhibits complex and 
highly nontrivial anomalous multiafhne dynamic scaling 
[]6|-pl| is an additional theoretical incentive in under- 
standing its growth universality. Another significant fea- 
ture is that, by construction, this growth model is the low 
temperature version of the full temperature depen- 
dent activated diffusion MBE growth model |4|j|] because 
in the limited mobility model only the adatoms without 
any lateral bonding are allowed to increase their coor- 
dination through diffusion, and therefore it has consid- 
erable experimental significance. It may be appropriate 
in this context to point out that several experimental 
measurements |^ of MBE growth exponents (/3 w 0.2, 
a ss 0.7) are consistent with the (d=2-|-l dimensional) 
critical exponents given by the fourth order nonlinear 
conserved growth equation which, as we show in this pa- 
per and in ref. 12, defines the universality class of the 
limited mobility discrete growth model of ref.l. 

This work has been supported by the US-ONR and the 
NSF-MRSEC. 
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